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Abstract. We present the deep UBVRI observations of 
the HDF-S NICMOS field obtained as part of the Sci- 
ence Verification of the VLT Unit 1 telescope. The im- 
ages have been used to construct object catalogs and to 
obtain photometric redshifts. The effective field of view 
is ~ 70 x 70 arcsec 2 , and the formal 5a limiting magni- 
tudes (in a 2 FWHM aperture) are 26.3, 27.8, 27.5, 26.9, 
25.2 in the U, B, V, R and I bands, respectively. Thanks 
to the sub-arcsecond image quality, relatively long expo- 
sure time, and large collecting area of the VLT, this is 
the deepest set of multicolor images ever obtained from a 
ground-based telescope. 

Galaxy counts have been derived independently in 
each band, and show no significant departures from pre- 
vious data from wider areas. A multicolor photometric 
catalog of all the galaxies selected in the R band has also 
been obtained and used to derive photometric redshifts 
for all galaxies with R < 26.5, using also the J,H and K 
magnitudes from the NICMOS deep observations. A sig- 
nificant fraction (~28%)of the galaxies is placed at z > 2. 
Among them, six robust candidates are found at z > 2.5 . 
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1. Introduction 

Deep imaging of extragalactic fields has long been recog- 
nized to be a powerful tool to understand galaxy evolution 
(see Ellis 1997 for an extensive review). Although faint 
galaxy counts have been of paramount importance to show 
that galaxies do evolve with redshift, the overall scenario 
and the physical processes that led to galaxy evolution are 
still debated. A new approach developed in recent years 
uses deep multicolor surveys to study the fainter magni- 
tude galaxies: deep multi-band images are taken with a 
comolete set of broad-band filters, in order to cover the 



populations at different redshifts. The Hubble Deep Field 
North (Williams et al. 1996) is the best-known examples 
of this kind of observations, but ground-based images have 
also been used, mainly to define sharp color criteria that 
select high-redshift galaxy candidates (Steidel et al 1995, 
Giallongo et al. 1998, Arnouts et al 1998). 

This paper analyses deep observations of the Hubble 
Deep Field-South (HDF-S, Williams et al 1999) obtained 
in five colors (U, B, V, R and I) in August 1998 as part 
of the Science Verification phase of the first VLT 8.2m 
telescope (UT1). A description of the Science Verification 
programme at the VLT is to be found at 
fittp: //www. eso.org/ paranal/sv/ 



The data have been taken for a field centered on 
a = 22/i32m51.7, 6 = -60d38'48.2" (J2000), thus 
providing the optical complement to the near-IR J, H 
and K band images obtained with the HST NICMOS 

instrument (Fruchtcr et al 1999, see also 

tittp: //www. stsci.edu/ftp/science/hdfsouth/hdfs. htm | ), 
since the area covered by the optical observations is 
somewhat more extended than the NICMOS field of view. 

This paper is organized as follows: In Section 2 we 
describe the data reduction procedures, while in Section 3 
we discuss how the photometric catalogs were constructed. 
In Section 4 we derive the photometric redshifts for all the 
objects brighter than R ~ 26.5 and discuss briefly their 
distribution at high redshift. 

2. The Data Sample 

2.1. Observations and data reduction 

The data used here were retrieved from the ESO pub- 
lic archive, and refer to the observations that were ob- 
tained in the period August 17-September 1, 1998, us- 
ing the VLT Test Camera (VLTTC) at the Cassegrain fo- 
cus of the UT1. The VLTTC is a simple imaging camera 
which reimages the focal plane onto a 2048 2 , 24 fj,m pixels, 
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Table 1. Summary of the Observational Data 



Filter 


Number of 


total 


Final PSF (FWHM) 


Field of View 


Zero-point 


Limiting Mag. 




frames 


exp. time (s) 


(arcsec) 


(arcsec) x (arcsec) 


(i) 


(5 a) « 


U 


16 


17800 


0.73 


68.8 x 72.8 


31.58 


26.3 


B 


15 


10200 


0.78 


71.7x71.7 


33.13 


27.8 


V 


16 


14400 


0.82 


78x69.5 


33.85 


27.5 


H 


8 


7200 


0.64 


72.6x 66.9 


34.06 


26.9 


I 


12 


10160 


0.71 


75.6 x 67.0 


33.44 


25.2 



(1) : Referred to the images released by ESO 

(2) : Magnitude limits are computed from the standard deviation of the sky counts inside apertures with diameters of 2x 
FWHM without aperture correction. 



program a set of several exposures through the standard 
U BVRI Johnson-Cousin filters was obtained with single 
exposure integration times ranging from 600 to 1200 sec- 
onds. Airmass of individual exposures was always < 1.4, 
with a median value of 1.25. Observations were obtained 
following the standard criteria for deep imaging, i.e. ap- 
plying a slight offset between individual pointings to allow 
the removal of the detector imprints. Full details on the 
Test Camera, the CCD detector, the filter curves and on 
the individual exposures are available on the WEB site 
http: / /www.hq.eso.org/paranal/sv. 

The data reduction has been carried out with different 
softwares at ESO and at the Rome Observatory. Although 
the two pipelines used completely independent environ- 
ments and tools, they are quite similar in concept. Most 
of the steps are identical to those described elsewhere (e.g. 
Giallongo et al. 1998; Arnouts et al. 1998), and will not be 
repeated here. It is worth noting that flat-fielding is par- 
ticularly critical in the VLTTC, since its CCD suffers from 
a very large blemish near the center and other lesser de- 
fects in the whole area. The central blemish is wavelength 
dependent, while variations from night to night caused by 
moving dust grains are also noticeable in the flats. Nev- 
ertheless, a satisfactory solution was finally found by con- 
structing a separate "super-flat" as the median image of 
the unregistered images for each observing night. The fi- 
nal accuracy in flat-fielding is estimated to be better than 
1%. 

Only frames with seeing better than 1 arcsec were used 
in the final coaddition, without applying any drizzling al- 
gorithm (Fruchter and Hook 1998) since the VLTTC sam- 
pling is always much smaller than the seeing. As is cus- 
tomary in dithered multiple exposures, the edges of the 
final images are of poorer quality, since only a limited 
number of frames contribute to the observed flux. In the 
present observations the problem is noticeable because a 
large dithering pattern had to be applied to remove the 
blemishes and because of the small size of the field of view. 
Moreover, the central pointing of the coadded frames was 
not exactly the same in the different bands, which has fur- 



regions independently, keeping only the inner regions that 
were covered with 100% of the observations. These images 
were used to extract independent catalogs in each band. 

We also produced a set of coadded images that are 
trimmed and aligned to the central field of the R frame, 
and these images were used to prepare the multicolor cat- 
alog. In practice, about 20% of this field is not covered 
by all the U and B frames, and therefore the coadded im- 
ages are slightly shallower in these bands. The recently 
available images obtained by NICMOS with the F110W, 
F160W and F222M filters (Fruchter et al 1999) were also 
rebinned and aligned to the VLTTC R frame. 

Table | summarizes the observational data that have 
been used, the FWHM of the PSF of the coadded single 
color images, their area, the photometric zero point, to- 
gether with the formal 5c limiting magnitudes. These were 
conservatively computed by taking the cr-clipped standard 
deviation of the sky counts in an aperture 2 x FWHM wide, 
taken at random positions on the images. 

2.2. Photometric calibration 

The photometric calibrations were obtained by reducing 
a series of standard stars from the Landolt (1992) sample. 
Self-consistent photometric solutions to the standard 
Johnson-Cousin system were derived for several individ- 
ual nights, along with average solutions that use all the 
data from photometric nights. The zero points, color 
terms and extinction coefficients derived are listed at 
fittp: //www, hq.eso.org/ paranal/sv/html/data/photom.txt . 
We emphasize however that the coadded images presented 
here cannot be calibrated directly using the average coef- 
ficients, since they are the average of different exposures 
taken under various conditions. 

A few isolated relatively bright stars have been selected 
in the field. Then, images obtained during each nights for 
which a photometric solution exists have been taken, and 
accurate magnitudes of the selected stars have been mea- 
sured in each of these images. The resulting instrumental 
magnitudes have been converted into Johnson magnitudes 
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Fig. 1. Differential galaxy number 
counts for the U, B, R and I bands. 
VLT data are shown with open 
(raw) and filled (corrected) squares. 
Error bars take into account Poisson 
error only. Other data taken from 
the literature are shown. U Band: 
empty circles from Hogg et al 1995; 
crosses from HDF-N (Williams et al 
1996). A correction Uj = U A b - 0.8 
has been applied to the latter data. 
B Band: (open) filled circles are 
(un-)corrected data from Metcalfe 
et al 1995a. Crosses from Arnouts 
et al 1998. R Band: Filled circles 
are raw Keck data from Smail et 
al 1995. / Band: Filled circles are 
from the NTT Deep Field (Arnouts 
et al 1998), crosses from HDF-N 
(Williams et al 1996). 



The same photometry has been applied to the summed 
images, and the final zero point has been chosen in order 
to reproduce the magnitudes of the selected stars. We es- 
timate that the final accuracy of this procedure is ~ 0.05 
mags in each band. Finally, the zero points were corrected 
for galactic absorption with E(B — V) — 0.02 (Burstein 
& Heiles 1982) with SU = 0.095, SB = 0.084, SV = 0.063, 
SR = 0.05, SI = 0.038. 

3. The Photometric Catalog 

The analysis of the two sets of images (the individual coad- 
ded images and those trimmed to the R band image) was 
performed using the SExtractor image analysis package 
(Bertin & Arnouts, 1996) . 

3.1. Galaxy Counts 

Within SExtractor, images were smoothed with a gaussian 
filter matching the seeing, and the detection threshold was 
chosen at 3cr of the background intensity in a contiguous 
area of 1 FWHM. Following Djorgovski et al. (1995), for 
each object both isophotal and aperture magnitudes (in a 



ture one. For fainter objects, an aperture correction to 5" 
has been estimated on bright stars and applied to cor- 
rect the 2 FWHM aperture magnitude. This procedure is 
strictly valid for star-like objects only, but has been shown 
to be a good approximation on deep images (Smail et al. 
1995). 

Bright stars have been excluded from the catalogs us- 
ing the CLASS_STAR parameter provided by SExtractor. 
A threshold CLASS_STAR< 0.9 has been set, on the basis 
of the comparison between ground-based and HST data 
(Arnouts et al 1998). At fainter magnitudes (e.g. R > 24) 
the neural network classifier does not work properly any- 
more, but stars are not expected to dominate the counts 
and therefore have been ignored. All the single-color cat- 
alogs are available on the web site 
[http://www.mporzio.astro.it/HIGrTZ . 

A correction for incompleteness has been estimated in 
each band as in Arnouts et al. (1998), accounting for both 
false detection and the non-detection of real objects. We 
warn the reader that the correction for incompleteness has 
to be taken with some caution, because of the very lim- 
ited size of the sample used as a reference. The raw and 
corrected counts in each band are shown in Fig ITJ, and 
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Table 2. Galaxy candidates at z > 2.5 



OBJ# 


x w 




R 


•^phot 


VLT 73 


629.4 


463.9 


26.13 


2.5 


VLT 116 


478.1 


237.6 


25.32 


2.5 


VLT 98 


244.4 


332.0 


25.44 


3.05 


VLT 27 


747.6 


619.7 


24.29 


3.45 


VLT 96 


722.2 


339.4 


25.07 


3.9* 


VLT 24 


64.7 


630.5 


25.64 


4.65* 



(1) : Referred to the R image released by ESO 
(*) : unresolved in the VLT images. 



telescope, thanks to the sub-arcsecond image quality, rel- 
atively long exposure time and the large collecting area of 
the VLT. 



3.2. The Multicolor Catalog 

The multicolor catalog has been obtained from the set of 
aligned images, taking the R frame as reference. Object 
detection and the measurement of the R magnitude have 
been performed on the R frame exactly as described in 
the previous section. Then colors have been measured in 
a fixed circular aperture of 14 pixels (corresponding to 2 
FWHM of the R frame), keeping the object position found 
on the R frame. To allow for the seeing difference among 
the coadded images, colors have been measured on images 
degraded to the 0.82" seeing of the V band image. The 
R < 26.5 subsample used for the following analysis con- 
sists of 91 galaxies and is available on the WEB 



http: / /www. mporzio.astro.it/HIGHZ. The R = 26.5 thresh 



old has been set in order to ensure good photometric accu- 
racy and meaningful colors on the whole sample. Obvious 
stars have been excluded using the CLASS_STAR param- 
eter (see above) in the R and - when available - in the 
NICMOS images. 

Fig ^ shows the color distribution of the R < 26.5 
sample in the V — I vs U — V plane. Overplotted on the 
observed colors are the evolutionary tracks of few galaxy 
templates as a function of redshift. They have been chosen 
to broadly encompass the most common spectral types, 
and are based on the synthetic models of Bruzual and 
Chariot (GISSEL library, 1996). It is clearly seen that a 
significant fraction of the faint galaxies has colors typical 
of star- forming galaxies at z ~ 2, as expected in these 
very deep images (Metcalfe et al. 1995b). Though most 
of these objects are probably too faint for a spectroscopic 
confirmation with FORS or ISAAC at the VLT, their na- 
ture and redshifts can be investigated further by means of 
a photometric redshift analysis. 

4. The Photometric Redshifts of the R < 26.5 
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Fig. 2. The color distribution of galaxies in the VLT field. 
Points refer to galaxies in the R < 26.5 sample. Lower lim- 
its are shown at the 3tr level. Circles marks the objects at 
z > 2.5 selected from the photometric redshift. Overplot- 
ted are the evolutionary tracks of three galaxy templates, 
computed with a Miller-Scalo IMF and the following pa- 
rameters: r = 0.3 Gyr, Zf orm = 7 and E(B-V)=0.1 (solid 
line); r = 5 Gyr, z forrn = 5 and E(B - V) = 0.2(dotted 
line); r — oo, Zf orrn — 5, no dust (dashed line). Points at 
z = 0.5, 1, 1.5, 2, 2.5, 3 are marked by a cross. 
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Fig. 4. The broad band energy dis- 
tributions of high redshift galaxy 
candidates. Object VLT 154 was 
found by Treu et al 1998. Other ob- 
jects are galaxies at z > 3 in the 
VLT field. Magnitudes have been 
converted into the AB system by 
Uab = U + .71, Bab = B — 0.09, 
Vab = V - 0.003, Rab = R + -18, 
Iab — I+A2. The best-fitting spec- 
trum is also shown. 



using a technique extensively described elsewhere (Gial- 
longo et al. 1998). In brief, we have computed the expected 
galaxy colors as a function of redshift for synthetic models 
in the GISSEL library, for an extensive variety of combina- 
tions of age, metallicity, IMF, and e- folding star formation 
time-scale. The reddening produced by internal dust in 
star-forming galaxies has then been added using the SMC 
extinction law (Pei 1992), along with the absorption pro- 
duced by hydrogen in the intergalactic medium (Madau 
1995). Finally, at any redshift galaxies are allowed to have 
any age smaller than the Hubble time at that redshift 
(Q = l and H — 50 km s _1 Mpc _1 have been adopted 
throughout the paper). 

The resulting large dataset includes ~5x 10 5 "simu- 
lated galaxies" , and a classical ^-minimization procedure 
has been applied to find the best-fitting spectral template 
to the observed colors. 

This procedure has been tested on 108 galaxies 
with spectroscopically confirmed redshifts in the HDF-N 
(Cowie 1997; Cohen et al. 1997; Dickinson et al. 1997; 
Lowenthal et al. 1997 Fernadez-Soto et al 1998), obtain- 
ing an accuracy a z ~ 0.1 in the redshift interval z = — 3.5 
(Giallongo et al. 1998; Fontana et al. 1999). 

The resulting redshift distribution is shown in Fig [| A 



A very red object was already identified 
candidate from the preliminary NICMOS observations, 
from being detected in the H band (F160W) and unde- 
tected in CTIO 4m telescope R and I band images with 
(R - H)ab > 3.9 and (I - H) AB > 3.5 (Treu et al 1998). 
This object (named VLT 154 in our catalog) is detected 
in all the VLT images except the U band, and our pho- 
tometric redshift analysis indicates a redshift z ~ 1.8 (see 
Fig ||) . The redshift accuracy is limited by the lack of ma- 
jor spectral features, and acceptable solutions in the range 
1.5 < z < 2.1 can be found using different combinations 
of the parameters involved, consistent with the result of 
Stiavelli et al. (1998). The z = 1.8 best-fit is obtained 
with solar metallicity, E(B - V) = 0.15 (with a SMC ex- 
tinction law), star-formation timescale r = 0.3 Gyr and 
an age of 2 Gyrs. Adopting the Calzetti (1997) extinction 
law we obtain z — 1.65 with E(B — V) — 0.2. Assuming 
no reddening from dust, we obtain an even higher best- 
fit redshift of z = 2.05. At z = 1.8 , this object would 
have Mk = —25.67 and Mb = —21.84 (fc-corrections 
are computed exactly from the best-fitting spectrum). 
Spectroscopic follow-up with ISAAC will hopefully re- 
veal whether this object is indeed a high redshift ellip- 
tical galaxy undergoing passive evolution, as suggested by 
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The list of the objects at z > 2.5 is given in Table g, 
while Fig. || shows the best fitting spectra of five of them. 
Two of these objects are unresolved in the VLT images, 
although they are too faint for the CLASS-STAR param- 
eter to be reliable. Since they also fall outside the HST- 
STIS image overlapping the HDF-S NICMOS field, we are 
not able to exclude the possibility that they are actually 
galactic stars, that are the major source of interlopers in 
the Steidel et al. (1996) sample. 

Before comparing these results to the HDF-N it is 
worth noting that HDF F300W filter is significantly wider 
and bluer than the Johnson U used here, with only a small 
overlapping region. As a result, the redshift range sam- 
pled by the HDF is wider and centered at a lower redshift 
(2 < z < 3.2) than in the present work (2.5 < z < 3.4), 
and the surface density of "U -dropout" galaxies in HDF- 
N with V606 < 26.5 is significantly higher than found here, 
or ~ 18 arcmin ~ 2 (Pozzetti et al. 1998). 

Only one galaxy candidate at z > 4 results from the 
redshift distribution in Fig. ||. This is consistent with the 
number density of 0.8 arcmin -2 "B-dropout" galaxies de- 
tected in the HDF-N (Pozzetti et al. 1998) and lower than 
a similar prediction in the NTT-Deep Field (2.7 arcmin -2 
at r <26, Arnouts et al 1998). 

Given the small size of the field studied with the VLT- 
TC these results are obviously of limited statistical signif- 
icance. They demonstrate however the potential of future, 
wide field VLT instruments like FORS-1 and FORS-2, now 
planned to become operational in 1999 and 2000. These in- 
struments will allow to explore to an unprecedented depth 
the distribution and evolutionary status of galaxies in the 
early universe. 
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